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Bose- Einstein condensates (BEC) have recently been the subject of considerable study as possible 
analogue models of general relativity. In particular it was shown that the propagation of phase 
perturbations in a BEC can, under certain conditions, closely mimic the dynamics of scalar quantum 
fields in curved spacetimes. In two previous articles gr-qc/0110036 , gr-qc/0305061 we noted that a 
varying scattering length in the BEC corresponds to a varying speed of light in the "effective metric" . 
^ | Recent experiments have indeed achieved a controlled tuning of the scattering length in Rubidium 

. 85. In this article we shall discuss the prospects for the use of this particular experimental effect 

to test some of the predictions of semiclassical quantum gravity, for instance, particle production 
in an expanding universe. We stress that these effects are generally much larger than the Hawking 
radiation expected from causal horizons, and so there are much better chances for their detection 
in the near future. 

PACS numbers: 04.70.Dy, 03.75.Fi, 04.80.-y, cond-mat/0307491 

> ■ 

I. INTRODUCTION — MOTIVATION 

^ ■ 

. Semiclassical gravity has played a central role in theoretical physics. Phenomena such as the Hawking effect or 
cosmological particle production are commonly considered to be crucial first steps on the way to building up a 
consistent fully-quantum theory of gravity (see for example pj). However a fundamental limit to these investigations 
is imposed by the fact that their most basic description is based on linear QFT on a fixed (classical) continuum 
spacetime. Several theoretical approaches have been developed to overcome this limitation: In a fashion that we can 
call "top-down" , string models [brane models] have in some special situations developed a high energy description of 
the Hawking effect 0, while "bottom-up" approaches, based on Stochastic gravity and the Einstein-Langevin analysis 
of particle creation by a gravitational field, have in recent years provided further insig ht@,i. 

On the other hand, the physics community has so far lacked any possibility for direct experimental tests of these 
ideas. Indeed this lack of experimental guidance is a severe hindrance with respect to further developments in 
O ■ semiclassical gravity (or full-fledged quantum gravity for that matter). In particular we have no experimental guidance 
regarding the manner in which the predictions of curved spacetime quantum field theory are changed once the 
hypotheses of non-discreteness and/or a non- fluctuating background are relaxed. In this regard the analogue models 
of gravity developed in recent years can be considered a first attempt to create an arena which can serve as a theoretical, 
jj^ • and possibly observational, laboratory to test aspects of these scenarios. 

No experimental set up has yet been realized in which the predictions of analogue models can be observationally 
tested. Nevertheless theoretical analyses of analogue models 0,0 have been so far remarkably successful in teaching us 
how semiclassical gravity phenomena are sensitive to possible quantum gravity effects, such as, for example, modified 
(Lorentz violating)dispersion relations p|. (See for example the trans-Planckian problem in the Hawking effect Q 
and in cosmology (g.) 

What we intend to discuss in this article is a particular class of experiment — that we hope could be realized in the 
very near future — wherein certain analogue gravity model predictions could be tested. The interest in doing so would 
not just be that of confirming a now well-established theoretical prediction, but mainly trying to evince deviations 
from the naive theoretical predictions due to the intrinsic discrete nature of the experimental system and/or to the 
possible role of non-linearities. 
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We shall focus our attention on the analogue gravity system established by the propagation of linearized phase 
perturbations in a Bose- Einstein condensate floL llll fill HH. IbH ] . In particular we shall consider an experiment 
where a time- varying scattering length is used to simulate the cosmological expansion of the universe, and its associated 
quantum creation of particles. 

It is interesting to note that in reference the authors proposed an explanation of the so called "bosenova" 
phenomenon 0] (a controlled instability of the bulk condensate induced by a sign variation of the scattering length) 
through a particular implementation of a version of analog cosmological particle production. In that approach the 
entire bulk of the condensate is rendered unstable and suffers catastrophic breakdown. Our current paper takes 
a complementary approach: Instead of trying to explain an observed phenomenon via analog cosmological particle 
production, we shall instead consider the most favourable conditions to observe it — preferably without violent 
disruption of the entire condensate. 

The scheme of the paper will be as follows: In the next Section we will review the physics of BECs regarding its 
analogue gravity features. Section III will be devoted to the discussion of how to simulate a FRW universe within 
a BEC. There exist two main routes to do this. One considers an explosive expansion of the condensate; the other 
makes use of the possibility of tuning the strength of the interaction between the different bosons in the condensate. 
This latter route will be the main concern of this paper. 

In Section IV we will first qualitatively describe how the modification of the interaction strength (encoded in the 
value of the scattering length) yields cosmological particle creation. Next, in Subsection A we will discuss [in the 
context of current BEC technology] whether there exists a regime in which this particle creation process can actually 
be reproduced. We will see that there is a limit on the rapidity of change of the background configuration, associated 
with the need to enforce a "Markov approximation" , in order for the whole construction to make sense. However, 
this bound still leaves a lot of parameter space available to look for the particle creation effect. Subsection B reviews 
the cosmological particle creation process, emphasizing the particular features of BEC systems. Then, Subsection C 
discuss the actual observability of the effect. Finally, we conclude with a short summary and discussion. 



II. ANALOGUE GRAVITY IN BOSE— EINSTEIN CONDENSATES 



Bose-Einstein condensates (BEC) have recently become subject of extensive study as possible analogue models of 
general relativity (3, lift 111! Il2t Il3t 1 1 11] . In particular it was shown that the propagation of phase perturbations in a 
BEC can under certain conditions closely mimic the dynamics of quantum fields in curved spacetimes. In previous 
papers we noted that a varying scattering length in the BEC system corresponds to a varying speed of light in the 
"effective metric" 0, ^| . Recent experiments have indeed achieved a controlled tuning of the scattering length in 
Rubidium 85 ^tJ- The effect is powerful enough to lead to large non-perturbative changes in the effective metric. Let 
us start by very briefly reviewing the derivation of the acoustic metric for a BEC system. 

In the dilute gas approximation, one can describe a Bose gas through a quantum field '5 satisfying 

ih 4* = (-^-^ 2 + V oxt (x) + K (a)&$) $. (1) 



dt \ 2m 

Here k parameterizes the strength of the interactions between the different bosons in the gas. It can be re-expressed 
in terms of the scattering length as 

. , Airah 2 

k a = . (2) 

m 

As usual, the quantum field can be separated into a macroscopic (classical) condensate and a fluctuation: = ip + p, 
with (\&) = ip. Then, by adopting the self-consistent mean field approximation [Tsj 

^ (pp ~ 2 (p^ p) (p + (p(p) p\ (3) 

one can arrive at the set of coupled equations: 

d ( h 2 \ 

ih —ip(t, x) = - — V 2 + Vext(x) + Kn c ) ip(t, x) + k (2hip(t, x) + mip*(t, x)) ; (4) 
at \ 2m / 



■n 9 - 



-— V + Vext(x) + K 2n T \(p + k m T <f 
2m 



(5) 
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Here 



n c = |^(*> x )l' 
riT = n c + n\ 



m c = i/j 2 {t,x); 
fa = (00); 
rriT = m c + m. 



(0) 
{') 
(8) 



The equation for the classical wave function of the condensate is closed only when the back-reaction effect due to 
the fluctuations are neglected. (This back-reaction is hiding in the parameters m and h.) This is the approximation 
contemplated by the Gross-Pitaevskii equation. In general one will have to solve both equations simultaneously. 
Adopting the Madelung representation for the wave function of the condensate 



V>(£,x) = y/n c (t,x.) exp[-i9(t,x)/h], 



(9) 



and defining an irrotational "velocity field" by v e V#/m, the Gross-Pitaevskii equation can be rewritten as a 
continuity equation plus an Euler equation: 



d 

—n c + V • (n c v) = 0, 
at 



m— v + V 

at 



n 2 v 2 ^r c 



2m 



= 0. 



(10) 

(11) 



These equations are completely equivalent to those of an irrotational and inviscid fluid apart from the existence of 
the so-called quantum potential V^ uan tum = — fi 2 V 2 y / 7i^/(2m v /ro^), which has the dimensions of an energy. Note that 



v jV^ uan tum — V^ 



h 2 V- 



2m 



= V, 



n c ViVi lnn c 

4m 



which justifies the introduction of the so-called quantum stress tensor 

h 2 



quantum 

4 



n r VjV,- lnn r . 
4m J 



(12) 



(13) 



This tensor has the dimensions of pressure, and may be viewed as an intrinsically quantum anisotropic pressure 
contributing to the Euler equation. If we write the mass density of the Madelung fluid as p = m n Cl and use the fact 
that the flow is irrotational then the Euler equation takes the form 



-v+(v.V)v 



P V 



V ext (t,x) 



up 



2m 2 



^ i ^quantum q 



(14) 



Note that the term V ex t/m has the dimensions of specific enthalpy, while up 2 /{2m) represents a bulk pressure. When 
the gradients in the density of the condensate are small one can neglect the quantum stress term leading to the standard 
hydrodynamic approximation. Because the flow is irrotational, the Euler equation is often more conveniently written 
in Hamilton-Jacobi form: 



d n 

m—6 + 
at 



(V6 



2m 



+ V ext (t,x.) + nn c 



h 2 v 1 



2m 



0. 



(15) 



Apart from the wave function of the condensate itself, we also have to account for the [typically small] quantum 
perturbations of the system (JSJ). These quantum perturbations can be described in several different ways, here we are 
interested in the "quantum acoustic representation" 



-w/h 



1 



n\ — i 



2^/rT c h 

where ni,9i are real quantum fields. By using this representation equation JSJ can be rewritten as 

d t ni + —V • (ni V6< + n c VoA = 0, 

m \ / 

d t 6i + —V6 ■ V0i + K(a) m - — D 2 n x = 0. 
m 2m 



(16) 

(17) 
(18) 



4 



Here D 2 represents a second-order differential operator obtained from linearizing the quantum potential. Explicitly: 



V 3/2 [V 2 (nt 1/2 )1 ni +\n- 



D 2 n x = --n- 3/2 [V 2 (n 



1/2 y2 (r 



-1/2 



ni). 



(19) 



2'"C L ' V"C /J '"J- 1 2" C * 

The equations we have just written can be obtained easily by linearizing the Gross-Pitaevskii equation around a 
classical solution: n c — > n c + n\ ,</>—*■</> + 4>i ■ It is important to realize that in those equations the back-reaction of 
the quantum fluctuations on the background solution has been assumed negligible. 

We also see in those equations, (|17ll and Q18[l. that time variations of V ex t and time variations of the scattering length 
a appear to act in very different ways. Whereas the external potential only influences the background equation lfT5|) 
[and hence the acoustic metric in the analog description] , the scattering length directly influences both the perturbation 
and background equations. 

From the previous equations for the linearized perturbations it is possible to derive a wave equation for 9i (or 
alternatively, for ni). All we need is to substitute in equation (|17|l the rii obtained from equation l|18|) . This leads 
to a PDE that is second-order in time derivatives but infinite order in space derivatives — to simplify things we can 
construct the symmetric 4x4 matrix 



f 



Hi 



(20) 



(Greek indices run from 0-3, while Roman indices run from 1-3.) Then, introducing (3 + 1) -dimensional space-time 
coordinates — x M = (t; x l ) — the wave equation for 9\ is easily rewritten as 

d^r &A) = o. 



Where the f^ v are differential operators acting on space only: 



r = 




k(cl) - 


h 2 

2m 


D 2 


-1 




fj = 




k(o) - 


h 2 

2m 




1 V J '6» 
m 


fiO = 


m 


k(o) 


h 2 ' 

- — D 2 

2m 


-l 



5 %i V^o 



k(o) 



2m 



Do 



(21) 

(22) 
(23) 
(24) 
(25) 



Now, if we make an spectral decomposition of the field B\ we can see that for wavelengths larger than h/mc s (this 
corresponds to the "healing length" , as we will explain below), the terms coming from the linearization of the quantum 
potential (the D 2 ) can be neglected in the previous expressions, in which case the f^ u can be approximated by numbers, 
instead of differential operators. (This is the heart of the acoustic approximation.) Then, by identifying 



-9 9 



{lis 



f 



I lis 



(26) 



the equation for the field Q\ becomes that of a (massless minimally coupled) quantum scalar field over a curved 
background 



A0i 



1 



d» {V^g f v d u ) 0i = o, 



(27) 



with an effective metric of the form 



m c s (a,n c 



-{c s (a,n c ) 2 - v 2 ) 



Here the magnitude c s (n Cy a) represents the speed of the phonons in the medium: 

n(a) n c 



c s (a,n c ) 



m 



(28) 



(29) 
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III. ANALOG MODELS FOR COSMOLOGICAL SPACETIMES 

To find analog models for cosmological spacetimes we will consider a generalized GP equation where the external 
potential and the coupling constant can both change with time 

= (-^-V 2 + V cxt (t,x) + K (t) |V(t,x)| 2 ^(i,x). (30) 

The technical steps in the calculation change in a straightforward manner and lead to a simple time-dependent acoustic 
metric 

ds 2 = [_( c 2 _ w 2) dt 2 _ 2 y.dxdt + dx 2 ] (31) 

it is this time-dependent effective metric that we now wish to use for simulating a cosmological spacetime and, 
subsequently probing cosmological particle production. 

A. Cosmological analog by explosion 

Starting from the geometry 131|) there are different ways in which one can reproduce a cosmologically expanding 
geometry. Following [pH Il9l I20L l2l| one can take a radial profile for the velocity v = (b/b)r, with b a scale factor 
depending only on t, and define a new radial coordinate as r& = r/b. In these new coordinates, the metric will be 
expressed as 

dg2 = _nc_ r2 dt 2 + b 2 {dr 2 + 2^2)1 _ (32) 

m c s 

Now, one solution for the BEC wave function that reproduces a FRW universe is this: Introducing a Hubble-like 
parameter, 

H b (t) = (33) 

the equation of continuity can be written as 

h c + 3H b (t) n c = 0; n c {t) = (34) 

with no ^constant. Then, the solution for 9 can be obtained from 115(1 

b 2 , f jj Kn o 



and it requires an external potential of the form 

W*,r) = -^ + ^mr 2 . (36) 

One could certainly construct such a potential in a "sufficiently large" region around r = 0; this would correspond to 
a "sufficiently large" part of a homogeneous and flat FRW universe. Outside this region, the potential will have in 
practice some confining walls. The final effective metric can be written as 

ds 2 = -T 2 (t) dt 2 + a 2 s (t) {dr 2 b + r 2 dn 2 ), (37) 

with 

,V4 . nl 1 



a s (t)=(^-b) , and T(t) ^ ^ / (38) 

V 772 K / k z at (t) 



(*)' 

So finally we end up with a FRW universe whose proper Friedmann time, r, is related to the laboratory time, t, by 
r = fT(t)dt. 
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Although this explosion route seems promising, one should note that this analog models has substantial drawbacks. 
In particular it is easy to see from equation (|35|l that one would get for the condensate a linearly rising velocity field 
v = V0/m oc r. Hence this particular realization of a FRW effective geometry is guaranteed to possess an apparent 
horizon, a spherical surface in which the speed of the fluid surpasses the speed of sound. From a dynamical point of 
view, this might introduce many practical problems not intrinsically inherent to the type of geometries we are trying 
to reproduce. Because of this, we view the use of an exploding medium as not being the preferred route for building 
an analogue for an expanding FRW universe. (For an alternate view, where the explosion route is the primary focus 
of attention, see 0.) 



B. Cosmological analog by varying speed of sound 



Another way to reproduce cosmologically expanding configurations in, wc think, a cleaner fashion is by taking 
advantage of the possibility to change the value of the scattering length offered by some BECs 01 • Let us note that 
in reference pH , the authors already used the existence of time dependence on the scattering length, in combination 
with a time dependent external potential, to create explosive configuration of the type described in previous subsection. 
Here, we are going to use this variability with a different strategy. 

Let us again start from H31|l but now with v — at all times: 

dg2 = _ncc 1 ^ + d _. 2 
m m c s 

Then, it is not difficult to envisage a situation in which n c is constant in a "sufficiently large" region (think of a 
"sufficiently large" close-to-hard- walled box). Then, the continuity equation is directly satisfied and the Hamilton- 
Jacobi equation tells us that with a fixed external potential the phase function 9 will depend only on t adapting 
itself to the changes of n(t). Changing the scattering length with time directly causes changes in the value of the 
propagation speed c s . (That temporal changes in the velocity of sound can be interpreted as a cosmological expansion 
without invoking any sort of velocity in the medium has already been suggested in in the context of superfluid 
Helium.) We now define r = J[n c c s (t) jm^-^dt and write 

ds 2 = -dr 2 + a 2 s (r) dx 2 , (40) 

where 

f n c \ 1 12 _ f n c \ 1/4 _ / n e \ 1/4 1 _ A 

as{T > ~ \mc s (T) ) ~{^(a)) ~ Vi^2 ) ^TJI=aU4- (41) 

In this model an expansion corresponds to a decrease of the scattering length and vice versa. The Hubble function 
for this spacetime is 

1 / i / \ 1/2 • 

1 c' 1 I m \ c 



■ 2 c s 2 \ n c ) c*J 2 



^ = - = -o- = -o - (42) 



(The prime represents derivative with respect Friedmann time, the dot derivative with respect laboratory time). This 
is the model we will consider in the following discussion. 



IV. ANALOG COSMOLOGICAL PARTICLE CREATION 



Let us now present a qualitative explanation of how is that we can closely simulate cosmological particle creation 
within this model. We can start with a condensate in a stationary state described, in a sufficiently large volume, by a 
constant background density n c and a phase function linear in time, = Eq t. This is a solution of the GP equation. 
For this, one needs to have a potential that reproduces a large enough close-to-hard- walled box and that satisfies 
the condition Eq = — V ex t — nn c . Apart from this classical background, there will be some quantum fluctuations 
over it. At temperatures much below the critical temperature these quantum fluctuations are very small and can be 
described by the Bogoliubov equations. (These quantum fluctuations are present even at zero temperature owing to 
the so-called quantum depletion phenomenon, see for example [23f .) Let us consider that these quantum fluctuations 
are in their vacuum state. If one now decreases the value of the scattering length in a sufficiently slow manner (this 
issue will be discussed bellow) , all the individual bosons (this is only an approximate concept in an interacting theory) 
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will be affected in the same way. This means that the value of the background magnitudes — these arc the coherent 
magnitudes — will be slowly modified. The GP equation tells you that the density function n c will be kept fixed while 
the phase function will develop a non-trivial dependence with time. At the same time the value of the speed of sound 
will decrease. 

Now, apart from the background configuration, what happens with the additional quantum fluctuations? The 
equation satisfied by the quantum fluctuations is, in the acoustic approximation (that is, for long wavelengths) that 
of a massless minimally coupled scalar field over an expanding background and therefore, it will yield cosmological 
pair production of particles. 

An interesting point to notice is that varying the external potential — by this we mean changing with time the value 
of the external potential, V cxt (t), in the central homogeneous region — changes the background configuration in the 
same manner as varying the scattering length does. However, as the speed of sound does not depend on the external 
potential, its variations will not lead to cosmological particle creation. 

Now, what can we say about the observability of the particle production process? The standard technique to 
"see" phonons is tomographic imaging. One opens the trap and looks at the expansion of the condensate. Phonons 
correspond to distributions for the momentum of the atoms in the trap and different momenta correspond to different 
travel distances of the atoms after you switch of the trap. Taking snapshots of this evolution shows the density 
contrasts and then the original momentum distribution. 

If the wave function \& is split into the condensate wave function tp plus quantum excitations (p (in this situation, 
we mean atoms) then, the density you observe is 

= + (43) 
as {$) — (tp) = 0. Therefore, the observability of the effect will depend on the value of the ratio 

or more simply on the spatially integrated ratio 

C=^-J 'dx 3 (fi®. (45) 

If this quantity is of order unity (say 1/2 or 1/10) then phonons can probably be "seen"; if it is 1/100 then seeing 
phonons is technologically difficult. If the particle production process was so strong that the calculation of C resulted 
on values close to unity or higher, this would indicate that the Bogoliubov mean field approximation would have been 
violated, and the BEC itself disrupted. 

Let us now perform some explicit calculations of the particle production expected in realistic situations with present- 
day BECs in which the scattering length a is changed in time from some initial value to a different final value. We first 
have to know how quickly we can drive these temporal changes while still ensuring that the different approximations 
involved in the analysis remain valid. 



A. Varying a, validity of the GP equation 

The previous analysis shows that in order to consider particle creation driven by a time- varying scattering length we 
must be sure to work in a regime where the background is "instantaneously" reacting to the changes in a. Moreover 
the very derivation of the effective metric description is based on the GP equation which we then want to make sure 
holds at each instant of time. 

So we must first determine the upper bound on the rapidity of the change in the scattering length a which still 
permits the GP equation to hold. This will also give an upper bound for the frequencies of the quasi-particles that 
might be created (if t is the shortest timescale over which we can drive the system then v w 1/r is the largest frequency 
of the quasi-particle we can create). The validity of the GP in describing the Bosc-Einstcin condensate is related to 
the validity of several crucial assumptions which permit us to perform certain approximations on the fundamental 
multi-particle Hamiltonian description. The relevant approximations are generally stated to be the "mean-field" 
approximation and the dilute gas approximation. It is nevertheless important to note that in a dynamical situation 
a third approximation, which we can call "Markovian" approximation, is implicitly assumed. 

Let us review the meaning and implications of these approximations: The mean field approximation is based on the 
assumption that most of the atoms are in the condensate phase and that the influence of the non condensed fraction 
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can be neglected. This implies that significant creation of quasi-particles with excessive energies can be dangerous. 
In particular from the Bogoliubov dispersion relation |24| 



h - 2 



c = y C 2fc2 + ^_ fc2 j 5 (46) 

we can deduce that excitation of quasi-particles with wavelengths comparable to the healing length would led to free 
particle states (for k > 27r/£; £ = h/(m c s ); u> ~ h 2 k 2 /(2m)). This argument seems to imply that one should require 
the typical timescale for the change in a to be no shorter than the healing time. (Which is the analogue in this 
situation of the Planck time in quantum gravity.) 

The dilute gas approximation is instead related to the way the interaction potential is simplified in the GP equation. 
This approximation is valid if n|a| 3 -C 1, so we shall have to keep the amplitude of the changes in na 3 small in order 
to satisfy this bound. We wish to emphasise that the dilute gas approximation does not appear to be a crucial 
approximation for the analog gravity picture to hold. As long as the interaction term K,\ip(t, x)| 2 can be generalized to 
be some higher order (but still local) term 7r(|-!/'|) an analog gravity description is not precluded (see, for example ^3). 

Finally the Markovian approximation is related to the fact that in dynamical situations the two-body time-dependent 
scattering matrix can have a complicated form due to the "memory" of the system (see e.g. section IV-A of the paper 
by Kohler and Burnett [25]]). Basically in these situations the system is described by a GP-like equation where the 
interaction term includes a "delay term" described by an integration over different times. The necessary assumption 
in order to have a Markovian description of the dynamics (which together with the two previous approximations 
leads to the GP equation) is then that the timcscales on which external parameters are changing are longer that the 
two-body collisional duration. That is, longer than the timescale over which a single interaction happens. (Reduced 
to the bare bones, we are asking that the scattering length does not change significantly during the period when a 
pair of atoms are interacting.) 

We can estimate the two-body collisional time by a simple calculation. All we need is the typical size of the region 
of strong interaction of two atoms in the condensate and the typical speed with which they move. The first quantity 
can be assumed to be of the order of the Van der Waals scale length: the inter-atomic potential V{r) is characterized 
by a short-range region of strong chemical bonding and a long-range Van der Waals potential, 

V -C 6 /r 6 . (47) 

This leads to a Van der Waals scale length 0, HHH [|J , 

Avdw = U^V /4 - (48) 



2\ h? 

This length is basically the size of the region of strong interaction: for r <C A v dw the scattering wave function oscillates 
rapidly due to the strong interaction potential. In alkali ground state interactions, Cq is the same for all hyperfine 
states of a given atomic pair; consequently, A vc iw is the same for all collision channels. For example in the case of 
Na2, it is about 2.4 nm. We shall assume here generically A v d\v ~ 1 nm. 

Regarding the typical speed of the atoms, this is set by the de Broglie momentum generated by the trap confinement: 
p = h/R and v = p/m. We shall assume a trap of typical size of ten microns. We then get 

A v dw _ A v dw m R 

tint — Z — j • 

v h 

We now confront this quantity with the timescale we have to be faster than in order to create modes with wavelengths 
shorter than the condensate size R. This is < S i ZC = R/c s . So 

^int A v dW m c s A v dW /rrA 

7 — = 7 = —r- ( 50 ) 

^size " S 

For typical BEC systems £ « 1/zm-O.l/xm (assuming that the scattering length ranges from one to a hundred nanome- 
ters) so 



t 



tint = AvdW Kl0 _3_ 10 - 2i (51) 



Note that t int can be computed to be 



VdW ^vdW 



m R 
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lCT 9 m • 85 • 10 9 eV • 1CT 5 m 
9-I0 16 m 2 /s 2 2tt • 6.5 • lO" 16 eV • s (52) 
w 2 ■ 1CT 6 s, 

so a microsecond is the shortest timescale allowed for the change in a. Note that this interaction time is shorter than 
the healing time ihcai = £/ c « ~ 10~ 3 — 10~ 5 s which plays the role of the "Planck time" in this system. Thus the 
GP equation is valid in the entire "sub-healing" regime, which is the primary regime of interest, and continues to 
hold well into the "trans-healing" regime (although the previous comments regarding the breakdown of the mean-field 
approximation in this regime remain valid). 



B. Analytical calculations: Changing a over a finite amount of time 



Now that we have estimated how fast the change in the scattering length can be driven, we can propose a particular 
time dependence and derive an estimate for the relative production of phonons. Particle production in an expanding 
universe has been extensively studied in the framework of semiclassical gravity (see e.g. |l|). In this regard the scope 
of this subsection will be to present an example of these calculations to non specialists as well as to evaluate the 
experimental feasibility of an experimental test. 

As a test-bed we shall choose a slightly simplified version of Parker's model [3(| • The FRW metric with flat spatial 
sections can be alternatively written as 



ds 2 = -a 6 s {r]) drj 2 + a 2 (77) dx 2 



(53) 



in which we are using for convenience a special type of pseudotime rj = J (n c /ma^)dt with a s — (ric/mcs) 1 ^ 2 . 

The scale factor is independent of x hence the mode decomposition for the quantum scalar field can be written as 



l(?7,x) 



dk 3 



(2tt) 3 / 2 



, a k e 



where the V'fe (v) are solutions of the equation 



d 2 ipk 
dr\ 2 



+ a 4 s (r)) k 2 ip k =0, 



that satisfy the normalization condition 



(54) 



(55) 



(56) 



Imagine now that the scale factor a s undergoes a hnite amount of expansion in a monotonic fashion. This means that 
the scale factor passes from and initial value a S i at r; = —00 to a hnal value a s f at r\ = +00 (or what is the same, that 
the scattering length passes from and initial value to a final value af). Spacetime can be approximated at early 
times and at late times respectively by two different Minkowski spacetimes. As is standard, we will assume that the 
quantum scalar field is initially in the vacuum state associated to the initial Minkowski spacetime. Then, we want to 
calculate the particle content of this state in the final Minkowski spacetime. 

There is a particularly convenient choice of function a s (r)) for which the physics is clear and for which one can 
exactly solve equation (|53|l : 



tanh 



(57) 



(This is a slight simplification of the model considered in Parker |3(j ; a variation of this model can be found in the 
book by Birrell and Davies |]| , see pp. 60 ff.) If we now impose as boundary condition for a solution that at 77 ~ —00 



Mv -> -00) = (2JW&) 1/2 expi-ika 2 ^]), 
one finds a particular set of exact solutions P, 

1 



(58) 



y/2k<? s 



: exp[ 



x 2-F1 ^1 + ikT] [a 2 f - a 2 , ; ]/2, +ikr] [a 2 s} - a 2 ,]/2; 1 - ikr] a 2 sl ; ~[1 + tanh(77/77 )]^ 



(59) 
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These solutions correspond to positive in-going modes. Similarly, the precise form of the exact solutions that at 
■q ~ +00 admit the asymptotic form 

tp k (r) -> +00) = (2ka 2 sf y 1 ' 2 exp(-ika 2 f r)), (60) 

(the positive out-going modes) is 



2ka 2 sf 



x 2-F1 ^1 - ikijo[a 2 sf - a 2 sl }/2, -ikr) [a 2 sf - a 2 sl }/2; 1 - ikrj a 2 sl ; i[l - tanh(?7/?7 )]^ 



(61) 



The Bogoliubov coefficients relating the early time (in-going) and late time (out-going) bases are then defined as 

ip™(ri 00) = [a k exp(~ika 2 sf r]) + (3 k exp(ika 2 sf i])}. (62) 

From this expression, and the exact solutions above, one can show [ll I30ll3ll l32| that 



a , = ^ju, st ^ /U s ' 1 ' K A ' s t' IRX) 

a^ + alri-ikMa^ + a^mri-ikMalf + al}^ [ > 

R _ 2a sf a sl T(-ikr) a 2 si ) T(+ikr l0 a 2 sf ) 



and that 



, _ sinh 2 (7rfc7 ?D [ag / + og i ]/2) 



l afc l = ■ U( U 2^ \f u 2"T' ( 65 ) 

sinn(7rfc77oaJJ smh(wkr]oa^) 



sinh (Trkr]o[a 2 sf - a 2 J/2) 
sinh(7r/c77 a^) smh^krjoa 2 .*) 



I^l 2 = . u r 2 , 2 v (66) 



These expressions are related to the Bogoliubov coefficients by [cf Birrell-Davies equation (3.93)] 

a(k in ,k out ) = a k S 3 (k in - k out ), (67) 



P(kim k ou t) — [3 k 5 (ki n "f~ kout)- 

The spectrum of particles in the final state is then 

diV 



d 3 fc D 



(3(k 

in 1 f^out , 



which gives 



dN 



\/3 k \ 2 S 3 (0) S^km + kout) d 3 k in = \(3 k \ 2 S-\0) 



(68) 



(69) 



(70) 



We now use the standard scattering theory result that a momentum-space delta function evaluated at zero is propor- 
tional to the volume of the "universe" , in this case the volume of the BEC, 



5(6) 



V 



(2tt) 3 



to sec that 



dN 

d. kn-ii. 



V 

(2*y 



■\(3 k \ 2 . 



(71) 



(72) 
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Equivalently 



dfc, 

and the total number of emitted phonons is 



<! V —lAI 2 k 2 out , (73) 



N =^l h\ 2 l-' 2 ^- (7-1! 



As a practical matter the integral will always be cut off at high momentum — most typically by the inverse timescale 
ri^ 1 over which the propagation speed changes, but if nothing else the integral cannot be trusted for momenta higher 
than that associated with the healing length fchoai = 2tt/£ for the reasons previously discussed (see Sec. lIV A|) . 

In order to gain a better understanding of the particle creation just described, it may be useful to study separately 
the two opposite regimes characterizing this phenomenon. In fact, for a given timescale of change, 770, driving the 
particle creation, one has a simplified description of the particle production when considering the case of modes with 
frequencies much smaller than 1 / 770 (sudden approximation) or much larger than that (adiabatic approximation). 
After a brief discussion of these regimes we shall deal with the full intermediate case. 



1. Sudden approximation 

A particularly useful approximation is to take the "sudden limit" . Mathematically this consists of taking a step 
function for the scale factor 

ai( V ) = 4 + 6(77) [aif-ait]. (75) 

Physically this means that one is considering that the change in a s (77) is driven more rapidly than the frequency band 
one is interested in. In this case, this means that the change is so fast that the entire acoustic regime is excited (the 
time rate is trans- healing), but sufficiently slowly that the GP equation still continues to hold (the time rate is still 
sub-interaction). However we shall still have to put in "by hand" a high momentum cut off, given by the healing 
length fc max = fchcai = 27r/£, because beyond this point we cannot trust the dispersion relation to remain on the 
acoustic branch implicit in our calculation. 

The relevant calculation can be preformed by simply considering the mathematical 770 — > limit in Parker's model. 
Indeed the a and j3 coefficients are now momentum independent and 

M = -£ -, (76) 

ZCL S f Q>si 



1°?/-°^ (77) 



As should be expected, particle production in this sudden limit depends only on the change in the scale factor. The 
particle production spectrum is now flat (more precisely, phase space limited) all the way up to the healing frequency. 
A rigorous result is that for any monotonic change in a s (rj) from a S i to a s t the magnitudes of the a and f3 coefficients 
are less than or equal to those calculated for the sudden approximation — consequently the sudden approximation 
provides an absolute upper bound on particle production. 
The number of phonons produced is 

N = 7^2 l/3| 2 *Lx = ^l/3| 2 (2V0 3 = y ^l/3| 2 . (78) 

That is 

The good news for current purposes is that this scales as (i?/£) 3 . Now the trap size R is of order 10 microns, while the 
healing length is in the range from 1 to 0.1 micrometers; thus (i?/£) 3 w 10 3 to 10 6 . A prefactor this big is desirable 
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in terms of producing an observable effect. As for a (scattering length), this can range from 100 nm to 1 nm, so the 
ratio di/df is up to w 100. Since a 2 (scale factor) is oc a -1 / 2 we have a 2 ^/a 2 si up to « 10. Therefore 



la l f 2 a } ? <2 (80) 



which is of order one — so there is no enormous suppression coming from the Bogoliubov factor. All in all, we estimate 
that N « 10 4 to 10 7 phonons can be produced in the sudden approximation. 



2. Adiabatic approximation 

In contrast, when the parameter r]oka 2 f is large compared with unity (that is, for large enough momenta), and 



provided a s f a s i, we have 



1 



exp(27r?7ofca 2 i ) 
so that the spectrum of phonons in the final state is 

dN 1 



dVd 3 k (2tt) 3 exp(27r?7ofca 2 2 ) - 1 ' 



(82) 



This is a correctly normalized Planckian [black body] distribution. We can associate a temperature T to the final 
phonon content produced by the expansion. Before doing that, let us write down some useful relations between the 
frequencies associated with the pseudotime 77 and with the laboratory time t. Asymptotically (either in the infinite 
past or in the infinite future) the relation between the times t and r\ is given by t ~ a 2 ry/c s + const ~ (ma^ / n c )rj + const 
so that 

■si k- (83) 

/ - c Bf k. (84) 

Here, the upper indices indicate with respect to which time a frequency is defined, and the sub-scripts i and / identify 
whether a particular magnitude is evaluated in the initial or in the final configuration. Now it is easy to see that the 
laboratory temperature that we would associate to the final configuration in the adiabatic approximation is 

T= 2^ hCsfV ° la7 ^ (85) 

In order to estimate this temperature we need to convert the time scale over which the scattering length changes, 
770 from the pseudo time to the laboratory time. To do this we define 



4 - 


n c 


■n 

UJ< = 




ma% 




n c 


77 


Csf 
a sf 


ma tf 



dt 

to = m ~r 

drj 



(86) 

r;=0 



For the particular profile, equation l|57|l . that associated with this model equation H86fl evaluates to 



m ( a% + at f \ ( a A + a 4 



If we use this relation between to and 770 , the temperature would be 

T= J_ n l a ji+<l . 
inks t a 2 sf a 2 si 

In this way we check the intuitive idea that in the expansion process one would create phonons with frequencies 
inversely related with the temporal scale of change of the configuration in laboratory. The number of phonons in the 
final state is 

it 2 {2v m a 2 sl Y %^nla% 64tt5 {c sf t Q f 1 ' 2 ! 
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Contrary to what happens in the sudden approximation, the adiabatic approximation is taking into account the 
rapidity with which the configuration changes. Thus, the total number of phonons calculated by trusting the adiabatic 
approximation throughout the whole range of frequencies is suppressed with respect to the sudden approximation 
calculation by a factor of 0.1 times (iheaiAo) 3 = 0.1 x [(healing time)/ (evolution time)] 3 — this factor consisting of a 
dimensionless number coming from the detailed expression for the integral, times the cubed ratio of the healing time 
in the condensate to the time over which the scattering length is forced to change. With to as 10 iheai we still get 
somewhere between 100 and 10 5 phonons. Once to w 100 th ea i this is reduced to somewhere between 10 _1 and 100 
phonons; but in this case we will also run into problems from finite volume effects — to is then comparable to the 
sound crossing time for the condensate and the momentum space delta functions appearing above are smeared out 
due to the finite volume of the condensate. (This point is carefully addressed in a rather different physical context 
in |33|. though many of the mathematical manipulations appearing therein are very similar to the present situation.) 

Note that as to — > theai the adiabatic approximation calculation still results in one order of magnitude less than the 
equivalent calculation with the sudden approximation. Moreover, in this case, the whole range of observable phonons 
(with frequencies between the healing frequency and the trap frequency) is beyond the strict range of applicability of 
the adiabatic approximation (remember that r)oka 2 j: 3> 1). Therefore, to be more precise one will have to make an 
intermediate analysis, in between the sudden and the adiabatic regimes. 

Let us estimate the value of the temperatures associated with the adiabatic approximation for temporal scales of 
change within the observable window. From equation l|88l) we can see that this temperature will be T ~ 10 _11 K • s-tj^ 1 , 
and so, for t between 10 -2 s (associated with the condensate size) and 10 -5 s (the shortest timescale compatible 
with the acoustic approximation associated with the initial configuration) will range from 1 nK to 1000 nK. From 
this estimate we can already see that, by modifying the scattering length on time scales close to the healing time, one 
could produce a bath of (almost) thermal phonons so energetic that even the mean field approximation might break 
down (causing the complete disruption of the condensate). 

When the peak frequency tends toward the healing frequency we see that the low frequency part of the observable 
spectrum will develop important departures from thermality. We shall now turn to the general intermediate case 
and describe the full spectrum of phonons created in our toy model. Then, we will discuss the observability of the 
cosmological particle creation effect in terms of the ratio C defined in equation . 



3. Intermediate regime 

Let us now consider an intermediate regime: Look at large momenta and ignore for the time being any cutoff arising 
from the interaction timescale or the healing length, then 

\(i\ 2 — > cxp[-47T77o k a 2 si ] as k — > oo, (90) 

so the particle spectrum is always exponentially suppressed at sufficiently high momenta. If we pick 770 to be longer 
than a healing time (so that we cannot use the sudden approximation) but still sufficiently short that we cannot use 
the adiabatic approximation, then we will need to retain the full form of |/3| 2 . The resulting spectrum of phonons is 

diV = 1 sinh 2 [TTr,ok(a 2 sf -a 2 sl )/2] 
dVA^k (27r) 3 sinh[7r?7ofca 2 y] sinh [^770 ka 2 sj \ ' 

and the total number of phonons produced is 

V Z" 00 o sinh 2 [71-770 fc (a 2 f - a 2 sl )/2] , s 

Zn* Jo smh[7r?7ofca^.J sinh [nrjo ka^ \ 

We can now consider the actual spectrum described by equation (|91|l by choosing plausible values for an experiment. 
However in order to get the spectrum that might be observed one has to convert the relevant quantities in equation 1911) 
to the laboratory counterparts. 

Using the expressions in equation (|84|l . and the relation between to and 770 given in equation l|87l) . and alternatively 
rewritten here as 

, K + <t \ ( <+«'.A 
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the number spectrum can be written as (we are reintroducing here the explicit in and out indices) 



dN 

d k n >,. 



V 



(2tt) s 



which implies 



sinh 



4 t 
, a sf UJ out 
nt — „4 j 



4 i 
a si UJ in 



II 



st 



sinh 



— 

a„, + a 



sinh 



2tt- 



* i 



- w £»*o 



<5 + &owt^ 



A- ^ 1 H 3 ? 



(94) 



d/V 



sinh 2 


n 4,4 

°& + a s 


/ 


Kut) 2 


2tt 2 

sinh 


a tf t 

a si "T" "s/ 


sinh 


r 2 9 

"Si "T" " S / 


C 3 



(95) 



Regarding the relative range of the scale factor a s , we have already seen that it can be deduced from the experimentally 
plausible range for the scattering length and we shall take a 2 f/a 2 si « 10. 

For the final value of the speed of sound c s / we shall take c s f ~ 10 mm/s. In fact also the range of the speed of 
sound can be determined from the scattering length. This can be reasonably be varied in the proximity of a Fcshbach 
resonance from 1 nm to 100 nm. For an experiment with heavy alkali atoms (e.g. rubidium) the speed of sound will 
then typically range from few mm/s to 10 mm/s. 



2e-05 



1.5e-05 



dX/dm 



1e-05 



5e-06 




Frcqucncy/KIIz 



FIG. 1: Number spectrum, equation 1951 . as phonons per 1000 cubic microns per unit out frequency. We have set c/ = lOmm/s 
and a 2 g j/a 2 si ~ 10. The typical timescale to is conservatively set equal to 10 _3 s. 

The total number of phonons emitted is given by 



TV = 



1 



V 



4 i 4 
a sf+ a si 



swh<[x(a 2 sf /a 2 si -l)/2} 



dx = 



1 V 



4 i 4 

atf + at 



16^ (c s/ t ) 3 V al f al ) J " sinn^^/a 2 ,)] sinh[x] "~ 16^ ( Cs/ * )3 { a 2 sf a 2 sl ) 1 sfl s * h 



(96) 

so qualitatively we have the same behaviour as in the adiabatic approximation, modulated by a dimensionless function 
of the ratio a s //a s j. (This expression will remain valid as long as 770 is longer than the healing time; at which stage 
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one should switch over to the sudden approximation.) Using hyperbolic trig identities the previous integral can be 
re- written as 



with 



F(z) 



1 



x 2 < coth(za;) coth(a;) — 1 



1 



sinh(zx) sinh(a;) 



dx 



F(1)=0; F(oo) = 



C(3) 



0.3005; and F(l/z) = z 3 F(z). 



(97) 



(98) 



For z > 1 the function F(z) quickly and smoothly approaches its asymptotic value. 

In closing our analysis of the intermediate regime we want now to compare the exact spectrum equation (19511 with 
the spectrum obtained within the adiabatic approximation. Equation (|82|) can be easily rewritten in the laboratory 
variables as 



V 



dN 



exp 



, a sf a si t + 



a si + a sf 



1 



The combined plot of and is shown in figure |3 



2e-05 
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1e-05 



5e-06 




(99) 



iTcqucncy/KIIz 

FIG. 2: Comparison of the exact number spectrum, equation l|95|l darker line in the graph, with the one obtained in the 
adiabatic approximation, equation 1991 . Both the spectra are shown as phonons per 1000 cubic microns per unit out frequency. 
In both the cases we have set c/ = lOmm/s and a 2 s f/a 2 si ~ 10. The typical timescale to is conservatively set equal to 10 _3 s. 



It is evident that apart from a minor discrepancy at lower frequencies the two plots basically coincide from the 
peak frequency (approximately 2 kHz) on. This is not so surprising given that the adiabatic approximation implies 



Tjoka 2 f 3> 1 which in laboratory variables corresponds to 



w o«t*o > F 



a 



(100) 



One can easily check that this inequality starts to be satisfied for frequencies of the order of a few kHz (and holds for 
any larger frequency). 
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C. Observability of the effect 

In this subsection we will calculate the relevant ratio C defined in equation Q45[l. We want C > 1/100 to make the 
phonons easily detectable. By using equation (|16|l we can deduce that 

C =W c J dx3 (i<^> + Y 2 ^ - ^^'^) ■ ( 101 ) 

In the acoustic regime and for the particular configurations we are looking at we know that Kn\ ~ —dt0\. Considering 
now that \ni(x,t),6i(y,t)] — ih 5 3 (x — y) or [9\(x, t),dt0\(y, t)] — ih k 5 3 '(x — y) — we can deduce this starting from 
the fundamental commutation relation [^(x, t), ^(y, t)} = 5 3 (x — y)) — we can arrive at 

c = ±J ^ + %m + \m) . (102) 

To calculate this average, we can expand the real field operator 0\ in terms of the creation and annihilation operators 
associated with the final configuration 

- I dk 'i^Mm^ °* '"****" + sa <103) 

and consider the appropriate particle content for the quantum state. If Nk is the average number of particles with 
momentum k in the quantum state, the previous magnitude can be expressed as 



c = 2^M J ^ N ^ + 2 ~W J JMk) Nkdk (104) 

1 I' mr 2 f 1 

' hw{k) N k dk 3 + 2 s - / T—rrrNkdk 3 . (105) 



2N C (m C 2) J w * N c J hiu(k) 

To reach this expression we have to substract the vacuum contribution. If we consider now the sudden approximation 
to calculate this rate (an upper limit to what one could get) we obtain 

C= +3-^-1^, (106) 



8 mc 2 s ftijJhc&i j N ( 

with TV the one in equation (J7SJ. Now fi^hcai — 10~ 10±1 eV. Instead mc 2 ~ atomic number x 10~ 13±1 eV ~ 10 _11±1 eV. 
The factor fiwheai/wc ~ 10, so we have that the relevant number C ~ AN/N C . But N c ~ 10 6 and N ~ 10 4 - 10 7 . 
This gives C ~ 4 x 10~ 2 — 40. This number is based on the sudden approximation and will be smaller in a more 
realistic calculation. However, remembering the discussion on the adiabatic approximation we know that, for temporal 
scales of change to of the order of the healing time, the actual coefficient C cannot be smaller than 0.1 times the 
previous estimate, i.e. C ~ 4 x 10~ 3 — 4. Therefore, by implementing in laboratory an expansive process with to in 
an intermediate regime, in between the healing times thcai = £/ 'c s associated with the initial configuration (to — 10 
seconds) and the final configuration (to — 10~ 3 seconds), one should be able to observe the effect. 



V. SUMMARY AND DISCUSSION 



In this work we have discussed the possibilities that BECs offer for simulating Lorcntzian geometries of the cos- 
mological type in the laboratory. There are two inequivalent paths one can follow in this task. The first one is 
based on provoking an expansive explosion in the condensate by changing with time the characteristic frequency of 
an isotropic and harmonic confining potential. This option implies that the velocity profile of the condensate acquires 
arbitrarily high values at large distances from the center. So there is always a sphere at which the velocity of the 
expanding BEC would surpass its sound velocity: A sonic horizon would be formed. In practice, due to the fact that 
physically realizable BECs are finite systems, one can only reproduce on them a portion of an expanding universe. 
Therefore one might argue that the sonic horizon would be formed outside the BEC. However, the velocity of sound 
in BECs is so small (a few mm/s) that in practical situations the sonic horizon will be formed well inside the system. 
Now the existence of sonic horizons is certainly interesting in its own right, but is not inherent to the simulation of 
cosmological spacetimes. Moreover, the plausible dynamical instabilities associated with their formation could mask 
the observation of purely cosmological effects. 
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The alternative path to simulating a cosmological geometry that we have pursued in this article is to take advantage 
of the possibility of varying the scattering length or, what is the same, the interaction strength between the atoms 
in the condensate. In this case, what we need is a confining potential with a sufficiently large almost flat minimum 
in which a portion of the condensate stays at rest. In this configuration, there is no formation of sonic horizons and 
thus we think it is (both conceptually and technically) a much cleaner path to follow. 

The description of the condensation phenomenon naturally involves the separation of the system into a "classical" 
wave function (the condensate part) and quantum fluctuations. In the acoustic approximation we can think of these 
quantum fluctuations as phonons over a classical background geometry, in this case, the analogue of a cosmological 
spacetime. Therefore, we can use the tuning of the scattering length to simulate not only a classical cosmologically 
expanding universe, but the quantum phenomenon of cosmological particle creation. We have analyzed this well 
known process by using a minor variant of Parker's model for a finite amount of expansion |3d| . Then, by working 
with numerical estimates appropriate to currently accessible BECs in dilute gases, we have presented an analysis of 
the feasibility of observing the effect in real experiments. 

We have seen that there is a more than plausible window for the observability of the effect with current technology. 
In current BECs the scattering length can easily be varied from a 100 nm to 1 nm. This produces an expansion in 
the geometrical scale factor of about three times. The temporal scale of change of the scattering length cannot be 
arbitrarily short. It has to be slower than the time scale in which the interaction between two atoms proceeds. We 
have calculated this time scale to be of the order of microseconds. However, we have also seen that, by driving the 
previous finite amount of expansion in a temporal scale of change of about fractions of millisecond, one could start 
to detect the presence of cosmological particle creation. From here one could shorten the time scale down to tens 
of microsecond progressively amplifying the expected effect. By the time one reaches tens of microseconds the effect 
would have been amplified by a factor of a hundred, (with timescale still above the interaction time), opening even 
the possibility of totally disrupting the condensate. 

The relevant temporal scales of change for the effect to be observable are of the order of the healing time in 
the condensate. Therefore we expect that apart from the phonon spectrum calculated by neglecting the modified 
dispersion relations at high energies, there will be also some production of quasiparticles. To observe the purely 
cosmological effect one would have to keep this quasiparticle production under a certain level; thus, the temporal 
scale of change should not be driven significantly beyond the healing time. 

In our analysis we have neglected finite volume effects. However we shall now show that these effects are insignificant 
for the typical condensate we considered here. The fractional change in the number of particles produced due to finite 
volume effects is expected to be of order l/(Kh ca iR) = £/(27ri?)=(cutoff wavelength) /(2nx size of the condensate). 
The ratio between the healing length and the BEC Thomas-Fermi radius can be expressed as a ratio between the 
harmonic trap length and the scattering length. 



For a harmonic oscillator length of about 10 microns, N c « 10 6 atoms and a scattering length of 1 nanometer one gets 
R w 17£. For a scattering length ten times larger (easily achievable with a Feshbach resonance and still compatible 
with N c a 3 <C 1) and « 1 /xm one would get R s» 100£. This implies that £ y /(2irR) < 1% and hence finite volume 
effects are negligible. 

To conclude, our analysis suggest that it should be already possible to observe the process of cosmological particle 
creation in BEC analogue systems, by changing the scattering length from an initial value of about 100 nm to a final 
value of about 1 nm on times scales shorter than milliseconds but larger than tens of microseconds. 
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